function diff = J_0_pp_fun(m,param)

    diff = -2.*param.c.*(1-param.alpha)./param.sigma^2./(param.gamma_R_2-param.gamma_R_1).* ...
        ((1./(1-param.alpha)-param.gamma_R_1).*m.^(param.gamma_R_1-2).*int_rep1_fun(m,param).*param.gamma_R_1.*(param.gamma_R_1-1)...
        +(1./(1-param.alpha)-param.gamma_R_1).*(param.zeta+delta_m(m,param))./m.^2.*param.gamma_R_1 ...
        -(1./(1-param.alpha)-param.gamma_R_2).*m.^(param.gamma_R_2-2).*int_rep2_fun(m,param).*param.gamma_R_2.*(param.gamma_R_2-1) ...    
        -(1./(1-param.alpha)-param.gamma_R_2).*(param.zeta+delta_m(m,param))./m.^2.*param.gamma_R_2 ...    
        -(param.gamma_R_2-param.gamma_R_1).*(param.zeta+delta_m(m,param))./m.^2 ...
        +(param.gamma_R_2-param.gamma_R_1).*delta_m_p(m,param)./m );

end